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We study the variation of the differential conductance G = dj/dV of a normal metal wire in 
a Superconductor/Normal metal heterostructure with a cross geometry under external microwave 
radiation applied to the superconducting parts. Our theoretical treatment is based on the quasiclas- 
sical Green's functions technique in the diffusive limit. Two limiting cases are considered: first, the 
limit of a weak proximity effect and low microwave frequency, second, the limit of a short dimension 
(short normal wire) and small irradiation amplitude. 
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I. INTRODUCTION 

Superconductor/normal metal (S/N) nanostructures, 
where the proximity effect (PE) plays an important role, 
have been studied very actively during last two decades. 
Interesting phenomena have been discovered in the 
course of these studies. Perhaps, the most remarkable 
one is an oscillatory dependence of the conductance of 
a normal wire attached to two superconductors which 
are incorporated into a superconducting loop^''^. This 
phenomenon was observed in the so-called "Andreev 
interferometers" , i.e. in multi-terminal SNS junctions 
(see^Ti^ as well as reviews2r~ and references therein). 
The reason for this oscillatory behavior of the differ- 
ential conductance G = dj/dV is a modification of the 
transport properties of the n wire due to the PE, i.e. due 
to the condensate induced in the n wire. The density of 
the induced condensate is very sensitive to an applied 
magnetic field H and oscillates with increasing H. 

Theorjiiir— was successful in explaining the exper- 
iments and predicting new phenomena, including the 
re-entrance of the conductance to the normal state in 
mesoscopic proximity conductors'^""— and transitions to 
the TT-state in the voltage-biased Andreev interferometers 
due to non-equilibrium effectsi^r— . The non-monotonic 
behavior of the conductance in SN point contacts 
and controllable nanostructrures has been observed in 
Refsj^""— , and the change of the sign of the critical 
Josephson current in multiterminal SNS junctions has 
been found in Refsi^i^. Many important results of the 
study of the SN mesoscopic structures are reviewed in 
Refs^-ia. 

The so-called 7r-states have also been realized in equi- 
librium Josephson SFS junctions with a ferromagnetic 
(F) layer between superconductors^^"^^ or in SIS Joseph- 
son junctions of high- Tc, d-wave superconductors^'*'^^. 

A number of new phenomena have been discovered in 



thin one-dimensional N and S wires^— (see also^^ for 
a recent review and references therein). 

Mesoscopic SNS structures proved to be a promising 
alternative to Superconducting Quantum Interference 
Devices (SQUIDs) for certain applications, including 
magnetic flux measurements and read-out of quantum 
bits (qubits)^ with a potential to achieve higher than 
state-of-the-art fidelity, sensitivity and read-out speed. 
To achieve such challenging aims extensive investigations 
of high frequency properties of S/N nanostructures on a 
scale similar to that of SQUIDs are in order. 

Studies undertaken to date concerned mainly the sta- 
tionary properties of S/N structures. Experimental data 
on S/N structures under microwave radiation appeared 
only recently^^i^. As to theoretical studies, one can 
mention two paper a^^'^^ where the ac impedance of a 
S/N structure was calculated. However, measuring the 
frequency dependence of the ac conductance is not an 
easy task. It is more convenient to measure a nonlinear 
dc response (dc conductance) to a microwave radiation. 
Recently, a numerical calculation of the dependence 
of the critical Josephson current Ic in SNS junction 
on the amplitude of an external ac radiation has been 
performed^. 

In this paper, using a simple model we calculate the 
dc conductance of a normal (n) wire in an S/N structure 
(cross geometry) as a function of the frequency and 
the amplitude of the external microwave radiation. We 
consider the limiting cases of a long and a short n wire 
and show that the response has a resonance peak at a 
frequency 57 close to Ss/h, where is the energy of a 
subgap in the n wire induced by the PE. Our theory 
predicts novel resonances and can help to optimize quan- 
tum devices based on hybrid SNS nanostructures^ii^. 

We employ the quasiclassical Green's function tech- 
nique in the diffusive limit. This means that we will 



2 



N 



-V 



n 



■L.r n 



N 



+V 



FIG. 1: Structure under consideration 



solve the Usadel equation'*^ for the retarded (advanced) 
Green's function ^^('^) and the corresponding equation 
for the Keldysh matrix function (section 2). First, 
a weak PE will be considered when the Usadel equation 
can be linearized (section 3). We calculate the dc con- 
ductance of the n wire in this limit, assuming that the 
frequency of the ac radiation is low (il ^ T). In sec- 
tion 4, the opposite limiting case of a short n wire will 
be analyzed for arbitrary frequencies Q. We present the 
frequency dependence of the correction to the dc conduc- 
tance caused by ac radiation. In section 5, we discuss the 
obtained results. 



II. MODEL AND BASIC EQUATIONS 

We consider an S/N structure shown in Fig.[T] It con- 
sists of a n wire or n film which connects two N and S 
reservoirs (n and N stand for a normal metal, S means a 
superconductor). The superconducting reservoirs may be 
connected by a superconducting contour. The transverse 
dimensions of the n wire are supposed to be smaller than 
characteristic dimensions of the problem, but larger than 
the Fermi wave length and the mean free path I (diffu- 
sive case). This implies that all quantities depend only 
on coordinates along the wire (the x— coordinate in the 
horizontal direction and the y— coordinate in the verti- 
cal direction). The dc voltage 2V is applied between the 
normal N reservoirs, and the phase difference 2if exists 
between the superconducting reservoirs. The phase ip is 
assumed to be time-dependent 



ip{t) ^ ipo + (pn cos{Q,t) 



(1) 



and related to the magnetic flux $ inside the supercon- 
ducting contour: (p{t) = 7r$(t)/$o with $(t) = H{t)S, 
where H{t) is an applied magnetic field and S is the area 
of the superconducting contour; that is, the magnetic 
field contains not only a constant component, but also 
an oscillating one. 



For simplicity, we assume the structure to be symmet- 
ric both in the horizontal and vertical directions. This 
implies, in particular, that the interface resistances RnN 
at X ^ zizLx 

are equal to each other (correspondingly, 
Rns{Ly) = Rns{—Ly)). Our aim is to calculate the dif- 
ferential dc conductance G between the N reservoirs 
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dj_ 

dV 



(2) 



as a function of the amplitude of the ac signal and 
the frequency Q. 

The calculations will be carried out on the basis of the 
well developed quasiclassical Green's functions technique 
(see the reviewsiSi^iri^) which successfully was applied to 
the study of S/N structures (see for example^"— i^iiii^il) . 
In this technique, all types of Green's functions (the "nor- 
mal" and Gor'kov's functions as well as the retarded, ad- 
vanced and Keldysh functions) are matrix elements of a 
4x4 matrix 



^'-l 9^ 



(3) 



where g^^"^^ are matrices of the retarded (advanced) 
Green's functions, and is a matrix of the Keldysh 
functions. The first matrices describe thermodynamical 
properties of the system (the density of states, super- 
current etc), whereas the matrix g^ is used to describe 
dissipative transport and nonequilibrium properties. 

The matrix q satisfies the normalization condition^ 



(gog) it,t')^S{t-t') 



(4) 



where "o" denotes the integral product {gog)(t,t') = 
J dti g{t,ti)-g{ti,t') and "■" is the conventional matrix 
product. The Fourier transform performed as g{e,e') = 



J dtdt' 



g{t,t') yields {gog) {£,£') = 2^5{£ - e') 



where now {gog) (e,e') = / ^ g(£, £i)-g(£i, e'). 

The matrix of Keldysh functions g^ can be expressed 
in terms of the matrices g^^^") and a matrix of distribu- 
tion functions F: 



gK ^ g%F - Fog^ 



(5) 



where the matrix F can be assumed to be diagonaU^: 

F = foF+ + f3F_. (6) 

Here tq is the identity matrix and the third Pauli 
matrix. The function F- describes the charge imbalance 
(premultiplied with the DOS and integrated over all en- 
ergies it gives the local voltage), while F^ characterizes 
the energy distribution of quasiparticles. 



Due to the general relation^ 

gA{e,e') = -f3g'''He\e)T3 



(7) 
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one can immediately calculate after finding the 
matrix g^. That means that knowing the matrices g^ 
and F we can determine all entries of g. 

The Green's functions in N and S reservoirs are as- 
sumed to have an equilibrium form corresponding to the 
voltages ±.V and phases ±(p{t). For example, the re- 
tarded (advanced) Green's functions in the upper S reser- 
voir are 



9s 



^^^\t,t') = s{tyg^^^\t-t'ys\t') 



(8) 



where S{t) = c:xjp[if3(p(t)/2\ is a unitary transformation 
matrix and the Fourier transform of 559^'' (i — t') equals 



e A 
-A -e 



(9) 



with = ±^(£±iO)2 

describes the BCS superconductor in the absence of 
phase. The retarded (advanced) Green's functions in the 
lower S reservoir are determined in the same way with 
the replacement tp{t) —ip{t). The matrix g^^^^ in the 
right (left) N reservoirs is equal to g^l^^'^ with A = 0, 

.R.(A) , . 

I.e. g^^ = ±T3. 

In the reservoirs the matrix F{t,t') can be repre- 
sented through the equilibrium distribution Feq{£) — 
tanh(e/2r) via Eq. ([8]) 



A^, i.e. the matrix 550"^'' 



F{t,t') = Sit)-F,,it~t')SHt') 



(10) 



The phase ip{t) in the upper S reservoir is given by 
Eq. ([1]), and for (pN{t) in the right TV reservoir, we have: 
'fiNit) = 2eVt. Therefore in the normal reservoir (at 
the right) the matrix distribution function has diagonal 
elements F/v(£)ii,22 = tanh (^^ (£ ± eV^)) , and can be 
written as F]y{e) — foF]\[+{e) + f^FN-is) 



FN±i£) = 2 



tanh ■ 



eV , , e-eV 

± tanli 

2T 



2T 



fll) 



Thus, all Green's functions in the reservoirs are defined 
above. 

Our task is to find the matrix g in the n wire. In the 
considered diffusive limit it obeys the equation^ 

+ ||-T3 + t{eVnit)g - geVnit')) - DV{goVg) = 0, 

(12) 

where 73 is a diagonal matrix with equal elements 
(73)11^22 = '^3, Vn is a local electrical potential in the 
n wire. We dropped the inelastic collision term suppos- 
ing that Eth = D/L"^^^ > t-"^^, where D is the diffusion 
coefficient, Lmax = max{La;^y} and Tinei is an inelastic 
scattering time. This equation is complemented by the 
boundary condition'*^ 



9°dx.vg\x.y=±L^,y = ±KN,s[g,gN,s]o 



(13) 



where kn,s = ^/C^crRnN.ns), RnN,nS are the nN and nS 
interface resistances per unit area and a is the conduc- 
tivity of the n wire. Here we introduced the commutator 
[9, 9n,s]o ^ go gN,s - 9n,s ° 9- The current in the n wire 
is determined by the formula 



j = i-Tr{f3-27r (5 0^:1,9)^2 (*'*)} 



(14) 



The matrix element (50 9^3)12 is the Keldysh compo- 



nent that equals {g o dxg)i2 



'0x9 



K 



9x9^ 



Even in a time-independent case, an analytical solu- 
tion of the problem can by found only under certain 
assumptions^""— i^liMTii^. In the considered case of a time- 
dependent phase difference, the problem becomes even 
more complicated. In order to solve the problem analyt- 
ically, we consider two limiting cases: a) weak proximity 
effect and slow phase variation in time; b) strong prox- 
imity effect in a short n wire and arbitrary frequency 
of the phase oscillations. 



III. WEAK PROXIMITY EFFECT; SLOW 
PHASE VARIATION 



In this section we will assume that the proximity ef- 
fect is weak and the phase difference (p{t) is almost con- 
stant in time. The latter assumption means that the 
frequency of phase variation satisfies the condition fl <C 
T/h. The weakness of the PE means that the anoma- 
lous (Gor'kov's) part f^'^^') of the retarded and advanced 
Green's functions in the n wire g^-^^^ = g^'^^^r^ + f^'^^^ 



can be assumed to be small 



|/«(^)|«1. 



(15) 



The matrix f^-^"^^ contains only off-diagonal elements. 
The diagonal part obtained from the normalization is 



9 



fs sa ± -r3 1 



1 



(16) 



Now we can linearize Eq. for the component 11(22), 
that is, for the retarded (advanced) Green's functions. 
Then we obtain a simple linear equation 



(17) 



where k^^^-* — yJ^2ie/D. The boundary conditions ([13 
for the matrices f^^-^^ acquire the form 



a./«(^)-2«:^./^(^) 



= 0; 



y=-Ly 



(18) 
(19) 
= 0; (20) 
= 0. (21) 
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As follows from Eq. ([8]) the functions gs, fs,ip are 

(22) 



{iT2 COS If + iTi sin If) A/ £_^'^^\ (23) 
We took into account that ip{t) is almost constant in time. 



One can show that the solution for in the horizontal 
part of the n wire is 

P = *T2/(X), (24) 

f{x) = C cosh{9xx/Lx) + sgn(a;) 5 sinh(6':ra;/i:r) (25) 

where sgn(a;) is the sign function. Dropping the index R 
of the quantities Kg, gsi^), the integration constants C 
and S can be written as 



'Ox cosh6'i, + rjv svnhOx) ■ rsAcoscp/S^^ 

' W) ' 

(26a) 

(rN cosh 6*3; + 0x sinh6la;) • rs A cos (p/^e 



V{e) 



(26b) 



where P(e) = {rsgs{£)dx + riqOy) cosh{0x + 6y) + 
(rsgsi^VN + 0x0y)smh{9x + Oy), vn^s = 2KN,sLx,y, 

Knowing the function f^{x), one can find the correc- 
tion to the conductance of the n wire due to the PE. 
In order to obtain the current, we take the (12) com- 
ponent (the Keldysh component) of Eq. p^ . multiply 
this component by fs and take the trace. In the Fourier 
representation we get (compare with Eq. (2) of Ref.*"') 



M^{e,ip,x)dxF^{x) =c{e,(p). 



(27) 



where the function M_(e, 93, a;) = 1 + i (fix) + f*{xyf 
determines the correction to the conductivity caused by 
the PE and c(e, ip) is an integration constant that is re- 
lated to the current: 



(7 f°° 



(28) 



It is determined from the boundary condition that can 
be obtained from Eq. ([T3]l 



M^{e,ip,x)dxF^{x) = c(e,^) = i^IFn- ~ F^{Lx)]. (29) 

where i'{e,ip) — U{1 + ^f{Lx)'^} is the density of states 
in the n wire near the nN interface. Finding F-{Lx) 
and c(e) from Eq. we obtain for the current density 
(compare with Eq. (13) of Rei^) 



de 



F, 



N- 



(30) 



Here Fn_ is defined according to Eq. PT|) . _R„ — Lx/cr is 
the resistance of the n wire of the length Lx in the normal 
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FIG. 2: Bias voltage dependence of the normalized variations 
of the resistance contributions SR„n /R„n and SRn/R„. Pa- 
rameter values: (p = 7r/3, Ly/Lx = 1, £jv/A = 2.5 • 10~^, 
es/A = 5 • 10-^ i?„/i?„iv = 1. 



state, and {M{e,(p)~'^) = {1/ Lx) J^" dx M_{e,ip,x)~'^ . 
The first term in the denominator is the nN interface 
resistance and the second term is the resistance of 
the (0, La;) section of the n wire modified by the PE. 
The expressions for the DOS v{e^ p) and the function 
(M(e, (^)^^) are given in appendix. 

For the differential conductance G = dj/dV at zero 
temperature we obtain 



GiV,ipit)) - {RnN/i^ + Rn{MieV,p)-^)}- 



(31) 



In Fig. [5] we show the dependence of the nN interface 
resistance variation 6RnN — RnN / ^ RnN and the resis- 
tance variation of the n wire 6Rn — Rn{M{eV, (p)~^)—Rn 
on the bias voltage V for a fixed phase difference. It 
can be seen that the 6RnN is either positive or negative 
depending on V, while SRn is always negative, i.e. the 
PE leads to voltage dependent changes of the interface 
resistance (caused by the changes of the DOS in the n 
wire) and to a decrease of the resistance of the n wire. 

The conductance variation SG — G{V,(p) — Gn, is 
shown in Fig. [3] for various values of RnN/RN, where 
Gn = {RnN + Rn}^^ is thc couductancc of the n-wire 
in the normal state. These results have been obtained 
earlierii*ii^>i&. 

We are interested in the dc conductance variation aver- 
aged in time: SGav = (r2/27r) J^''^^ dtSG{V,(p{t)). First, 
from Eqs. (1251126^ we can extract the dependence of the 
function / on the phase (p: f{x, p) = f{x, 0) cos p. Hence 
we obtain M-{e,p,x) — 1 + dM^{e,0,x) cos'^ p where 
SAI^{e,p,x) = M-{e,p,x) — 1. At the same time, 
^{e, p) — I + Si'{e, 0) cos^ p with 6i'{s, p) = v{£, pi) — 1. 
These observations lead to the relation 



5G{V, p{t)) = 6G{V, 0) cos^ p{t) 



(32) 
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FIG. 3: Bias voltage dependence of the normalized con- 
ductance variation 5G/Gn- Parameter values: tp = 7r/3, 
Ly/L^ = 1, ejv/A = 2.5 • 10"^ es/A = 5 • 10"^ Different 
cases: a) R„/R„n = 0.5, b) R„/R„n = 1, c) R„/R„n = 2. 



IV. STRONG PE; SHORT NORMAL WIRE 

In this section we analyze the limiting case of a short 
n wire when the Thouless energy Exh — D/L^ is much 
larger than characteristic energies: Eth ^ Dnj^ g, T, eV. 
In this case all the functions in Eq. (|12p are almost con- 
stant in space and we can integrate this equation from 
{a;, y} = ±0 to {a:, y} = iL^^y over x and y coordinates. 
The term f^-dtg + dfg-Ts (in the Fourier representation 
—ief^-g + ie'g-T^) is considered as a constant and the 
term with the voltage V is omitted because we neglect 
the voltage drop over the n wire; the voltage drops across 
the nN, nS interfaces and is set to be zero in the n wire. 
Performing this integration and summing up the results, 
we obtain 

2(i, + Ly)A = X{L,) - J,(-L,) - 4(+0) + 4(-0) 

+ Jy{Ly) - Jy{-Ly) - J,( + 0) + Jy{~0) (34) 




FIG. 4: Bias voltage dependence of the normalized time- 
averaged conductance variation SGav/Gn- Parameter values: 
ipo = 7r/3, Ly/L^ = 1, en/ A = 2.5 • 10'^, es/A = 5 • 10-^ 



eV/A = 5-W-\ R^/Rnt 



which by averaging over time yields 



SGa. = 5GiV, 0) • - (1 + Jo(2^o) cos(2(^o)) (33) 



where Jq is the Bessel function of the first kind and 
zeroth order. This oscillatory behavior of the time- 
averaged (dc) conductance variation SGav as a function 
of the ac amplitude can be seen in Fig. 21 

Thus, the calculations carried out in this section un- 
der assumption of adiabatic phase variations allow us to 
obtain the dependence of the conductance change SGav 
on the amplitude ipQ, but provide no information about 
the frequency dependence of SGav This dependence will 
be found in the next section. 



where Xix) ^ Dgod^gly^o, Jy{y) = D godyg\x=o, and 
A = —i{ef^-g — g-f^e'). Integration around the point 
(x,2/) = (0,0) yields the conservation of the "currents" 
(using terminology of the circuit theory^) 



4(+0) + Jy[+Q) = 4(-0) + Jy(-O) 



(35) 



Combining Eqs. p4ll35|) and the boundary condi- 
tions p^ . we arrive at the equation 



e'^s-ff - g-Tz£ = «£Af[5,.9Ar+]o + ies[g,gs-\ 



(36) 



Here eN,s = D/ (2RnN,nS<^L) is a characteristic energy 
related to the interface transparencies, L = + Ly. 
The energy e^r determines the damping in the spectrum 
of the n wire and the energy £5 is related to a subgap 
induced in the n wire due to the PE. The matrices gN,s+ 
are equal to: gN.s+ = ^[gNsiLx.y) + gN,si-Lx,y)]- 

In the limit of the short n wire considered in this 
section, we need to find only the retarded (advanced) 
Green's functions. Indeed, let us rewrite the expression 
for the current (|T4)) using the boundary condition (fT3)) 
at the right nN interface and concentrating on the dc 
component of the current: 



J = 



00 

— Tr{f3- / de ig'^g§+g''-gi-g§: 
IbeRnN J 



-9Ni 



where g^^"^^ 



(37) 

i^Fm^. The 



±f3 and Tr{f3-(g«-5^)} 
distribution function Fat- in the N reservoir is defined 
in Eq. (|lip . The integral over energies from the second 
and third term is zero because it is proportional to the 
voltage in the n wire which is set to be zero. Therefore 
the current can be written as 



1 



2eRnN 



de u{e)FM-{£) 



(38) 
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where = ^{g^ - g^) = 5R{g-'^(e)}. This formula 
has an obvious physical meaning - the current through 
the nN interface is determined by the product of the 
DOS in the n wire and N reservoir (i/^v = 1) and the 
distribution function in the N reservoir (the distribution 
function F- in the n wire is zero). 



Using Eqs. (g)), dTT]), (|3H1) we arrive at the following 
expression for the differential conductance: 



oo 



1 



1 



cosh 



2 e+eV 
2T 



cosh 



2 e-eV 
2T 



(39) 



In order to find the matrix g^, we can write the (11) 
component of Eq. (j36]) in the form 

(40) 



where e = e + ien, e' 



According to Eqs. (U]), ([S]) the matrix is a function 
of two times, gg_^_{t, t'), that is, in the Fourier representa- 
tion it is function of two energies: e,e'. Therefore, to find 
the matrix ^^(e, e') in a general case is a formidable task. 

However, we can assume that the amplitude of the ac 
component of the phase ipQ is small and obtain the solu- 
tion making an expansion in powers of ip^: 



9''^9E + 5w" + 62g^ + ^ 



(41) 



Here and later all matrix Green's functions written 
without arguments are functions of two energies (e,e')- 
Those of them which are diagonal in energy may be also 
(explicitly) written with a single energy argument, e.g. 

.9so+ =5so+(£)27r<5(e-£')- 



Similar to Eq. (1411) we represent the matrix (up to 

R 



the second order in iy9o) as g^^ = .950+ + + ^29s+ 



and find from Eq. ([SJ for the stationary part ggQ_^_ and the 
corrections Sig^^ (first order in (p^) and S2g§^ (second 
order in ipn): 



§§0+ = 27r(5o (efs + lA cos (^07^2) ^ 



(42) 



Si9§+ = -if2f</PfiAsin(^o(C ) {^^ + ^-^) (43) 

hfs+ = -*T2f A cos ^0 (Po + P2) (44) 

where we used the notation 8^ = (5(e — e' + w), = 
and defined the functions 

Po = <5o(2e-i + e,-;i^ + e-ij, 
P2 = \ (^20 + 8-2^) i^;,' + + (45) 

Using the expressions for Sig^^ and (5235+ given 
above we can calculate the corrections to gQ up to the 



second order in ipQ and the corresponding modification 
of the DOS 5v in the n wire. 

In the zcroth-order approximation, i.e. for ip^ = we 
obtain from Eq. (gO]) go {e,e') = g^[e)2-K5{e - e') where 
the matrix g^ie) obeys the equation 



[f:,Ef + iTAg^]^^ 



(46) 



containing Ef = e + i£s9soi^) = £ + + *£s5so(^)' 

The solution of this equation is^^ 

9^ie)^T,g^ie) + ^hf,"{e); 
g^{e) = E^Kf, f^ie)^EfUC^ (47) 



The quantity Ef^ 



where Cf = ^ [Ef f - {E^^f 
determines a subgap induced in the n wire due to the 
PE. Indeed, consider the most interesting case of small 
energies assuming that {£,£5} ■C A; then, « lA, 

/lo(^) ~ -* a'^d Cf ~ \/{e + ieNY - {sscostfoY. This 
means that the spectrum of the n wire has the same 
form as in the BCS superconductor with a damping 
£jv and a subgap £5 |cos(^oIj which depends on the nS 
interface transparency and phase difference. 

Note that the formula for the subgap induced in the 
N metal due to the PE in a tunnel SIN junction was 
obtained by McMillan^. 

The obtained results for the functions g^{e) and f^{s) 
can be easily generalized for the case of asymmetric nS 
interfaces with different interface resistances RnSi.2 (cor- 
respondingly, £51,2)- In the limit £si_2 ^ A, we obtain 
for the subgap Egg 



£sg{Vo) = \ \J £si + £|2 + 2£Sl£S2 COS 2ipQ 



(48) 



This formula shows that that the subgap as a function 
of the phase difference varies from 5|£S2 — £si| for 
tpo = tt/2 to \{es2 + £si) for tpo = 0. 

We proceed finding the corrections of the first {5ig^) 
and second {52g^) order in (^n for g^ in a way similar to 
the one used in^i^. The correction of the first order 5ig 
(for brevity, we drop the index R) obeys the equation 



Cego{£)-5ig - Sig-go{e')Ce' = i£s[9o,Si9s^ 



(49) 



which contains all terms of the first order in ipQ 
from Eq. (l40t . Note, that we are making use of the 
relation goie) — C,'^^(£T'i iesgsQ+ie)) evident from 
Eqs. (gni), (I4B1I471) . 



In order to solve Eq. (|49l) . it is useful to employ the 
normalization condition for g = ~g^ which for the 
first-order term oi g o g yields 



go{£)-6ig + 6ig-go{e') = 



(50) 



7 



From Eqs. (j49ll50|) . we find 

Si9s+ - 9o{£)-Sigs+-go{£') 



Ce + Ce 



(51) 



We determine the correction in the same manner. 
Reading off the second-order terms in Eq. PO]) gives 

[Cego, S2g]o = ies{[go. S2gs+]o + [Sig, Sigs+]o^ (52) 
The second-order part of the normalization condition is 
go{e)-62g + 52g-go{e') = -5ig o 5ig (53) 
Thus, we obtain the second-order correction 
hgs+ - ffo(£)-'525s+-5o(e') 



les- 



Ce + Ce' 

[(5i5s+,(5i.g]o Ce{^ig°5i9) 



Ce + Ce 



Ce + Ce 



■go{e'). (54) 



In order to calculate the correction to the dc conduc- 
tance caused by the ac radiation, 5G, we need to find 
Tr{f3-(5i5} and Tr{f3-(525} and take their parts propor- 
tional to 2TrS{e — e'). By inspection of Eqs. (j43|), (|5T|) one 
recognizes that the first-order correction contains only 
terms proportional to S{e — e' ± i7) and therefore only 
contributes to the ac current. This is the fundamental 



reason why the second-order analysis is needed to 
determine the variation of the dc conductance. 

As a result we just have to find the multiple of 
2T:d{e — e') contained in "Ti{T^-52g} which we denote as 
25dcg{£), that is (5dc.g(£)27r(5(e - e') := \ Trjfg-cSzgldc- 

We represent the function 6dcg{s) as a sum 



5dcg{s) = 5f^g{s) + 8^^^g{s) 



where the function S^i^ g{£) originates from the first term 

in Eq. ([54]) and the function S^^^ g{e) from the second 
and third terms. If we consider the case when the subgap 
Es |cos(/?o| is much less than A, i.e. 



(55) 



es |cosv5o| < A 



(56) 



then, at low energies e ^ es, the function g{e) is 
almost independent of fl, whereas the function S'^'^ g{e) 



depends strongly on at £ ss eg |cos(/3o| 
validity of Eq. (|56|) we obtain 



Assuming the 



-iel^^cos^o^, (57) 



J 



/ 1,2 2 • 2 .9o(£)[l + /o(£)/o(£ + ^) + .9o(£)go(£- 



[Ce + Ce+n? 



-m , .foie)[goie)foie + n) + foie)goie + n)] 



(58) 



r 



where the functions 50(e), /o(e), Ce are defined in 
Eq. (gT]). The sum sign index "±57" in Eq. dSS]) means 
that the given expression is added to the same one with 
the negative frequency {—^)- 

Using the function 5dcg{e, ^) we can calculate a cor- 
rection to the DOS 5v{e, il) due to the PE and with the 
aid of Eq. (HH) find the correction 5G{V, VL) to the differ- 
ential dc conductance. As follows from Eq. (P^ . at zero 
temperature the normalized differential dc conductance 
G{V, n) = G{V, n)RnN is equal to 

G{V,n) = Gq{V) + 5G{V,VL) = VQ{eV)+5v{eV,n) (59) 

with the definitions V{){eV) — K(go(eT^)) and 
5v{eV,Vl) = ^{5dcg{eV,n)). 

Using the obtained formula for (70(e) and 5dcg{,e) we 
can calculate the conductance Go and its correction 
5G due to microwave radiation for different values of 
parameters (damping e^r, phase difference 2lpq, etc.). 
The dependence of the conductance in the absence of 



radiation Go versus the applied voltage V is presented 
in Fig. O We see that this dependence follows the 
energy dependence of a SIN junction. In our case the n 
wire with an induced subgap plays a role of "S" with a 
damping £jv in the " superconductor" . Since the value 
of the induced subgap, e^g = £s \ cos 1^0 1 , depends on the 
phase difference 2lpq, the position of the peak in the 
dependence Go(V') is shifted downwards with increasing 
LPQ. Note that in an asymmetrical system (£51 ^ £52) 
the lowest value of the subgap is not zero (cf. Eq. 

In Fig. [5] we show the bias voltage dependence of 
the conductance correction due to ac radiation 5dcG 
(coefficient in front of ipo) for different values of (po- The 
magnitude and the position of the arising peaks depend 
strongly on the values of the parameters, e.g. Lp^. 

By varying the stationary phase difference tpQ or the 
damping e^ one can change the frequency dependence 
of the correction 5dcG considerably. This is shown in 
Fig. [7] and Fig. |5] respectively. One can see that if 
£jv ^ esg('/?o): then the dependence 5dcG(^) has a peak 
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Go 




eV 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 A 

FIG. 5: Normalized stationary differential conductance Go 
versus bias voltage V. Parameter values: T/A = 10~^, 
£s/A = 0.1, En/A = 10~^. Different cases: a) ~ tt/S, 
b) (fio = 7r/4, c) ifio = 37r/8. 

SG 



9l 




FIG. 6: Normalized second-order correction of differential 
conductance 5G versus bias voltage V . Parameter values: 
T/A = 10-^ es/A = 0.1, en/A = 10-^ Q./A = 5 ■ 10"^ 
Different cases: a) <^o = 7i"/8, b) tpo — n/'i, c) <^o = Svr/S. 

located at w esg((/3o) and split into two subpeaks. The 
splitting becomes more and more distinct with increasing 
bias voltage V. With decreasing ip^ and increasing etv, 
the form of this dependence changes significantly. For 
example, the resonance curve becomes broader with 
increasing damping. Increasing temperature leads to a 
similar effect as one can see in Fig. [9l 

In Fig. [TU] we plot the normalized conductance correc- 
tion 5dcG{ipo) as a function of for different values of 
the bias voltage V . At large V this dependence is close 
to sinusoidal, but at smaller voltages the form of the pe- 
riodic function 6dcG{ipQ) becomes more complicated. 



V. CONCLUSION 

We have calculated the change of the conductance 
in an S/N structure of the cross geometry under the 



SG 




0.00 0.05 0.10 0.15 0.20b. 



FIG. 7: Normalized second-order correction of differential 
conductance &G versus ac frequency fi. Parameter values: 
T/A = 2 • 10-^ £s/A = 0.1, EN I A = 10-^ el//A = 10"^ 
Different cases: a) ipo = tt/B, b) ipo ~ /'^, c) (/?o = 7r/3. 

6G 




0.00 0.05 0.10 0.15 0.20A 



FIG. 8: Normalized second-order correction of differential 
conductance SG versus ac frequency il. Parameter values: 
T/A = 2 ■ 10"^ es/A = 0.1, eV/A = 10"^ = 7r/3. 
Different cases: a) ejv/A = 5 • 10"^, b) en/A = 10"^, 
c) en/ A = 2 • 10"^ 

influence of microwave radiation. The calculations have 
been carried out on the basis of quasiclassical Green's 
functions in the diffusive limit. Two different limiting 
cases have been considered: a) a weak proximity effect 
and low frequency of radiation; b) a strong proximity 
effect and small amplitude of radiation. 

In the case a), the conductance change SG consists 
of two parts. One is related to a change of the nN 
interface resistance due to a modification of the DOS 
of the n wire. At small applied voltages Vjv, it is 
negative. Another part is caused by a modification of 
the conductance of the n wire due to the PE. This part 
is positive and consists of two competing contributions. 
One contribution, which is negative, stems from the 
a modification of the DOS of the n wire. Another 
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5G 




0.00 0.05 0.10 0.15 0.20 A 



We assumed that the nS interface resistance is 
larger than the resistance of the n wire, that is: 
pL <C RnS- This inequahty can be written in the form 
es <C ETh = D/L"^, that is, the subgap energy in the 
n wire is much smaller than the Thouless energy Sth- 
In the opposite limit, es 3> Srh, a gap of the order of 
ETh is induced in the n wire. This limit can be studied 
numerically. However, one can expect that in this limit 
the resonance should take place at uJres ~ £Th/fi- Exper- 
iment performed in Reff2£ corresponds to this limit. The 
frequency corresponding to the Thouless energy in ex- 
periment is equal to STh/h « ^(400/10~^)s~^ w 6GHz, 
whereas the resonance frequency is Vres ~ lOGHz. 



FIG. 9: Normalized second-order correction of differential 
conductance SG versus ac frequency Q, Parameter values: 
es/A = 0.1, en/A = 10-^ eV/A = 10-^ ipo = tt/3. Differ- 
ent cases: a) T = 0, b) T/A = 2 • 10"^ c) T/A = 6 • 10"^ 
d) T/A = 10'^ 
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FIG. 10: Normalized second-order correction of differential 
conductance 5G versus stationary phase difference ^po. Pa- 
rameter values: T/A = 2 • 10"^, es/A = 0.1, £n/A = 10"^ 
Q./A = 10"^ Different cases: a) V = 0, b) eV/A = 2 • 10"^ 
c) ey/A = 4-10-^ 



contribution is positive and caused by a term, which is 
similar to the Maki-Thompson termi^i^. The conduc- 
tance change 6G oscillates and decays with increasing 
amplitude of radiation. 



Note that we considered a simplified model. For exam- 
ple, we did not account for the change of the distribution 
function in the n film (heating effects). One can give 
estimations when the "heating" can be neglected. The 
change of an "effective" electron temperature 5T in the 
n wire is approximately given by: 5T « Te-phO'E'^ / Ce, 
where r^^ph is the electron-phonon inelastic scattering 
time, E < SYsRh/RbL = hn{RL/eLRb)fn is the 
ac electric field in the n wire and Ce ~ T ■ n/sp is 
the heat capacity of electron gas with concentration 
n. Taking into account that {Ri,a)^^ k, Z^//, we 
find that 5T/T < {eQ/T)^{Te^ph/T)Z^ipf^, where 
is the dimensionless coefficient of electron penetration 
through the SN interface, which is assumed to be small, 
I ~ VT is the mean free path in the n wire. There- 
fore, the heating would be very small if the condition 
< {eo/T)y/{T/re-ph)Z-^ is fulfilled. 

The obtained results are useful for understanding the 
response of the considered and analogous SN systems to 
microwave radiation which can be used, for example, in 
Q-bits. 

We would like to thank SFB 491 for financial support. 
One of us (VTP) was supported by the British EPSRC 
(Grant EP/F016891/1). 



In the case b) a short n wire was considered so that the 
resistance of the n wire is negligible in comparison with 
the resistance of the nN interface. The correction 6G 
has been found under assumption of a small amplitude 
of the radiation. We found that at small applied voltages 
V, the dependence SG{fl) has a resonance form. It has 
a maximum when the frequency Q is of the order of 
£sg = es I COS (fio I where Egg is a subgap in the spectrum 
of the n wire induced by the PE. With increasing V, the 
peak in the dependence SG{fl) splits into two peaks and 
overall form of this dependence becomes complicated. 



VI. APPENDIX 

The DOS in Eqs. (I29ll30p is given by the formula 
i^{e, if) = K[l + 5/(ix)^] with f{L^) defined in Eq. ([25]). 

Making use of the weak proximity effect approximation 
the function we rewrite {M{e, ip)~^) in Eq. (|30|) as 

(M(£,^)-i) = l-i(5?{/2 + /n) (60) 
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Using Eq. (|25l) one can easily calculate where 9'^ and 0" are the real and imaginary parts of 9^ 

respectively, i.e. 9^ — 9'^ + i9'J.. 

,,,, C^ + S^smh29,. , C^-S^ , ^^ sinh^d. 



{fn = 



2 2 

|Cp + |5|^sinh2g; \C\^ ~ \S\^ sm29'^ 
2 26i; ^ 2 29'^ 

,2 /)/ 2 /)// 



f / sinh 9' sin 0" \ 1 
< C*S' ( ^ + i - ) > (62) We use these expressions for calculating the function 

I V J ] (Af(£, (^)^^) and conductance G. 
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